Phase transitions in a network with range dependent connection probability 
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We consider a one-dimensional network in which the nodes at Euclidean distance / can have long 
range connections with a probabilty P(t) ~ l~ s in addition to nearest neighbour connections. This 
system has been shown to exhibit small world behaviour for 8 < 2 above which its behaviour is like 
a regular lattice. From the study of the clustering coefficients, we show that there is a transition to 
a random network at 5 = 1. The finite size scaling analysis of the clustering coefficients obtained 
from numerical simulations indicate that a continuous phase transition occurs at this point. Using 
these results, we find that the two transitions occurring in this network can be detected in any 
dimension by the behaviour of a single quantity, the average bond length. The phase transitions in 
all dimensions are non-trivial in nature. 
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A network is a collection of nodes which are connected 
either directly or indirectly by links. There are two ex- 
treme examples of networks: the regular and the ran- 
dom network. In a regular network (with infinitely many 
nodes), the probability that any two arbitrary nodes are 
connected is vanishingly small while in random networks, 
this probability remains finite. The two main properties 
which distinguish these two networks are the chemical 
distance and the clustering coefficient. The chemical dis- 
tance is the average shortest distance between any two 
nodes and is a long distance property. In a network 
with L nodes, the chemical distance typically behaves 
as Sl ~ hxL when it is random, while Si ~ L 1 ^ in a 
regular network in d dimensions. The clustering coeffi- 
cient C is the average fraction of connected triplets. Since 
the clustering coefficient measures the local connectivity 
structure it is a short range property. Typically C is high 
for the regular network and low for the random network. 

Recently, another kind of network, the small world net- 
work [Q , was proposed which shows random network-like 
properties at large scales and regular network-like prop- 
erties at small scales. Precisely, the chemical distance Sl 
behaves as ln(L) while C assumes a high value (compa- 
rable to a regular network) in this network. Small-world 
effect can be developed out of a regular lattice having 
local connections when long range links or connections 
are allowed to exist even with a very small probability. 

The underlying structure of a wide range of networks 
including social, biological, internet and transport net- 
works has been argued to be small world like 0. Addi- 
tional random long range connections in model systems 
like Ising chains or percolation networks also lead to new 
critical behaviour ^ ||, |j| . Many of these networks also 
show scale-free behaviour, i.e., if Q{k) is the number of 
nodes having k connections, then Q(k) ~ k f in a scale 
free network. 

In the Watts-Strogatz (WS) model 0, the nodes are 
arranged in a ring. Small world behaviour is observed 
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FIG. 1: A section of the network under consideration in 
(a) one and (b) two dimensions: the solid lines denote the 
nearest neighbours links between the nodes (circles). The 
dashed lines denote the long range connections occurring with 
the probability given in (hi). In (a), the nodes ABC form a 
triplet cluster. 



when nearest neighbour links are rewired randomly with 
a probability p. Later it was shown that there is a contin- 
uous phase transition occurring at p — > Q from a reg- 
ular to a small world phase. For all values of < p < 1, 
the network remains small world-like with Sl ~ hi(L) 
and a high value of C. Only at p — 1, the network be- 
haves like a random one when C vanishes. The transitions 
in the WS model, which is a standard prototype model 
for small world behaviour, are therefore trivial as the 
critical points do not separate phases of different critical 
behaviour. In a critical network with nontrivial phase 
transitions the critical points should have different phases 
on either side such that the small world phase emerges 
as a truly intermediate phase in between the random and 
regular phases. 
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In this communication, we show the existence of a 
network in which such non-trivial transitions can be 
achieved by tuning an appropriate parameter. In this 
network || |9|, [n^, in addition to the nearest neighbour 
links, connection to a node at an Euclidean distance I is 
present with a probability 



p(z)~r 



(i) 



Such a network can model the behaviour of linear poly- 
mers in which connections to further neighbours indeed 
occur with a power law probability. However, in order 
to simulate the polymer properly, additional constraints 
have to be considered . Studies on the physical layout 
of the Internet also strongly suggest that the connections 
are governed by a power law probability Jll|] . 

This network shows small world effect below 5 — 2 for 
the one dimensional case ||, 10 while according to ||, 
there is evidence of small world behaviour at 8 = d for 
the d-dimensional network. 

We have shown that there are two transitions occur- 
ring in such a network in one dimension: 

(i) A random to small world network transition at S = 1 

(ii) A small world to regular network transition at 8 = 2. 
The second transition is already known to exist ||, [l0| al- 
though the nature of the transition has not been explored 
in detail. We focus our attention on the first transition 
which has not been investigated before in either this or 
any other network to the best of our knowledge. We have 
been able to locate the critical point and the character- 
istic critical exponents for this transition. 

The other important result of our study is the iden- 
tification of a single quantity which can detect both the 
transitions. This analysis can be extended to higher di- 
mensions and the critical points are detected in any di- 
mension d. The phase diagram in the 6 — d plane shows 
that the phase transitions are non-trivial in all dimen- 
sions. 

First we examine the properties of this network in de- 
tail in one dimension. The one-dimensional network con- 
sists of L nodes occupying the sites of a one dimensional 
lattice of length L such that the nearest neighbour links 
are of unit length (see Fig. 1). In addition to the near- 
est neighbour links there are links between nodes at dis- 
tance I with the probability given by (|l|). Note that in 
one dimension, the Euclidean distance between two nodes 
coincides with the number of nearest neighbour links sep- 
arating them. The nearest neighbours links are always 
present in this network. This corresponds to the situa- 
tion in many realistic networks like linear polymers, Ising 
models etc. When the probability is normalised one has 



E 1< K i P(0 = 1. 



(2) 



The above normalisation condition essentially implies 
that the average number of long distance connections for 
each node is one. This enables a restriction on the net- 
work as the number of long distance bonds is conserved. 



If the value of 5 is made very high, most of the long dis- 
tance connections will be restricted to the near neigh- 
bours and one will effectively get a model with short 
range connections only. As 5 is made smaller, further 
neighbour bonds will be chosen. For very small values of 
S, connections to nodes at all possible distances are made 
and the network behaves like a random network. 

From the above picture and the knowledge of the exis- 
tence of a non-trivial phase transition from small world 
to regular network- like behaviour occurring at 6 — 2, we 
expect that all three kinds of behaviour, regular, ran- 
dom and small world will be present in this model - or 
in other words, there will be three regions along the 8 
axis: 5 < 8 C where it behaves like a random network, 

2) 
(1) 



8^ < 6 < 6 C 2 ^ where it is small world like and 6 > 8, 



where it is like a regular network. A non-zero value of 5c 
will signify a non-trivial transition. In order to study the 
transition from the random to the small world phase it is 
suffcient to study the clustering coefficients (the chemical 
distance has similar scaling behaviour in both phases). 
Here we consider clusters which are triplets with three 
members A, B and C. The condition that they form a 
cluster is that if B is connected to A and C, there is also 
a connection between A and C. Since the nearest neigh- 
bours are always present, we classify the possible clusters 
in three classes: 

1. Clusters with two nearest neighbour (nn) links each 
of length unity and one next nearest neighbour link of 
length 2. 

2. Clusters with one nn link, the other two links are of 
length li > 1 and l\ + 1. 

3. Clusters in which there is no nn link; the links have 
lengths h > 1, h > 1, and h + l 2 . 

Note that the triangular inequality of the link lengths 
is not valid in one dimension as the distances are always 
measured along the chain. 

Let Ci be the probability of the occurrence of a cluster 
belonging to the ith class (z = 1,2,3). In the continuum 
limit when I varies continuously, Ci take the forms: 



Ci 
C 2 

c 3 



JP(l)dl 

P{h)p{h + l)dh 

(fP(l)dl) 2 
J h h P(h)P(l 2 )P(h + l 2 )dhdl 2 

~" (fP(l)dl) 3 



(3a) 
(3b) 
(3c) 



In the last equation, the integration variables h, l 2 satisfy 
h +l 2 < L - 2. 

C = Sf =1 Ci is then the clustering coefficient of the 
system. We find that for 8 < 1, Ci vanishes for all i 
as L —> 00. In particular, to the leading order, C\ ~ 
L 5 - 1 while C 2 and C 3 are 0(L _1 ). The vanishing of C 
below 8=1 indicates that in this region the network is 
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Average bond lengths: I 
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Scaling plot for clustering coefficients 



0.9 r 

0.8 

0. 









6 £ = 2000 O 

- L = 1000 + 

'' L — 500 □ 

4 - L = 300 x 

3 _ L = 200 A 

2 - 
1 




[] 



-4 



-3 



-2 



1 

(6-l)L» 



FIG. 2: The data collapse of the scaled clustering coefficients 
are shown for one dimensional chains of different lengths L 
with periodic boundary conditions. Here a — 0.228 and b = 
0.258 



Average bond length: H 



3.4*x**0.21 
0.55*x**0.6 ■ 
0.24*x**0.78 ■ 



o o o 
A AAA 



random. For 8 > 1, all the three quantities remain finite 
in the same limit. Hence C may be interpreted as an 
order parameter like quantity which vanishes at 8 = 1 
where the transition from random to small world phase 
takes place. 

In order to find out the nature of the transition occur- 
ring at 8 = 1 (which we identify as 5^), we perform nu- 
merical simulations and find out the clustering coefficient 
for chains of different lengths with long range connections 
existing with a probability given by (|lj) . 



In thermodynamic system, the finite size scaling for a 
quantity $ is given by 

$(t) = L^"f(tL^ v ) (4) 

where t is the deviation from the critical point, <f> is the 
critical exponent associated with $ ($ ~ t -< ^) and v is 
the correlation length exponent. We use a finite size scal- 
ing form for the clustering coefficients which is analogous 
to that used in thermodynamic phase transitions: 

C = L-"g((SW - S )L b ). (5) 

Using 8 C = lj the value obtained in the continuum case, 
we find a very good data collapse when CL a is plotted 
against (6 - l)L b with the values a = 0.228 ± 0.017 and 
b = 0.258 ± 0.030 (see Fig. 2). These values are obtained 
using the Bhattacharjee-Seno method of data collapse 
|Q. For large values of x, it is expected that g(x) ~ x a ^ b . 
Note that the value of a/b is the estimate of the exponent 
/? as we have interpreted C as the "order parameter" (i.e., 
C ~ (5 - 8 { c ] ) p for L — > oo), therefore (3 = 0.89 ± 0.04. 
Also, the finite size scaling form indicates that v = l/b = 
3.87 ±0.51. 



FIG. 3: The scaling behaviour of the average bond length 
(l) for one dimensional chains of different lengths L is shown. 
In (a), (I) varies as L for S < 1. The curves are drawn for 
S = 0.2, 0.6 and 0.8 from top to bottom. In (b) the best 
fit lines for the curves for 1 < 8 < 2 show agreement with 
(/) ~ L 2 ~ s . The data are shown for 5 = 1.2, 1.4, 1.8, 2.2 and 
2.4 from top to bottom. For values of 8 > 2, (l) does not 
depend on L. 

From the above, we conclude that there is a continuous 
phase transition occurring at 6^ — 1 with characteristic 
exponents v ~ 3.87 and (3 = 0.89. Here (3 describes how 
the clustering coefficient vanishes as one approaches the 
random network and the exponent v is associated with a 
diverging length scale. Hence in this network the small 
world phase appears as an intermediate phase between 
the random and regular phases and the characterisric be- 
haviour of the network can be controlled by tuning the 
parameter 8. 



In general, the clustering coefficient and the chemical 
distances can detect either of the two phase transitions 
occurring in a network. This is because they have similar 
scaling behaviour in two of the phases and a different be- 
haviour in the third. Interestingly, in the present model, 
we find a quantity, the average bond length, which shows 
different scaling behaviour in each of these three phases. 
In the continuum limit, the average bond length (I) shows 
the following scaling behaviour: 

(I) - L {8<1) (6a) 

- L 2 ~ s (1< 8 < 2) (6b) 

- 0(1) (S>2). (6c) 

We immediately notice that the crossovers occur at 8^ 



4 



RANDOM 








/ sw. 








REGULAR 


1 2 


3 


4 5 



FIG. 4: The phase diagram of the network in the 8 — d plane. 
SW denotes the small world phase. 

(2) 

and 5c ■ Hence we find that this is a key quantity since 
both the transitions can be detected from it. This quan- 
tity is also simple to calculate. The numerical values of 
(I) for discrete lattices agree with the above results as 
shown in Fig. 3. 

Although the transition points can be located from the 
behaviour of (I), estimating the exponents is not straight- 
forward as it is difficult to cast the behaviour of (I) in a 
standard finite size scaling form as in (^). 

The one dimensional network which we have discussed 
so far can be generalised to any dimension d where the 
nodes occupy the sites of a d-dimensional hypercube of 
length L 1 ^ (see Fig lb). Each node is connected to its 
2d nearest neighbours and long range links to further 
neighbours occur with a probability given in (|l|). (Note 
that the structure of this network is different from that 
of a linear polymer embedded in a d dimensional lattice 
although the long range bonds occur with a power law 
probability in both systems.) For d > 1, neither 5c nor 

(2) 

S c is known. In principle, these can be estimated by 
studying the chemical distance and the clustering coeffi- 
cients ias in one dimension but the calculations become Electronic 



scale-free. On the other hand, when distance dependence 
in the form of ([l]) is introduced in a growing network, 
several interesting features are observed [fl3j . 

In summary, our analysis of a model network where 
the additional long range bonds are present with a prob- 
ability dependent on the Euclidean distance separating 
the nodes shows that there is a continuous phase transi- 
tion occurring at a finite value of the parameter 5 where 
the clustering coefficient behaves like an order parame- 
ter. The transition separates a random and a small world 
phase. With evidence of a transition from a small world 
to regular behaviour already existing, we find that this 
network can be tuned to show regular, random and small 
world behaviour for different values of 5. The different 
behaviour occurring on the two sides of the critical points 
mark the existence of non-trivial phase transitions. This 
is a feature absent in the familiar WS model of small 
world network, where the transitions are reminiscent of 
the zero temperature phase transitions occurring in the 
one dimensional Ising model. 

Comparing the results obtained from the clustering co- 
efficients and the chemical distances, we find that both 
the transitions can be detected from the behaviour of the 
average bond lengths. This analysis can be extended to 
any dimension and the transition points located. We be- 
lieve that this idea could be useful in general for locating 
critical points in network whenever the bond length is a 
meaningful quantity. 
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however, can be easily calculated in any general dimen- 
sion and we find (I) again shows different kinds of be- 
haviour in the three regions < 5 < d (where (I) ~ 
L),d < 5 < d+ 1 (where (I) ~ L d ~ s+1 ), and 5 > 
d+1 (where (I) ~ 0(1)). From the results of one di- 
mension we proceed to claim that the transition points 



in d dimensions are 5, 



(i) 



d and 5, 



(2) 



1. This 



indicates that the small world region again exists as an 
intermediate phase between a growing random region and 
a regular region. The width of the small world phase is 
independent of the dimensionality. The phase diagram 
in the 5 — d plane is shown in Fig. 4. 

As mentioned earlier, many real networks exhibit scale- 
free behaviour which implies that the degree distribution 
is a power law. Here we checked thai there is no scale free 
behaviour in any regime in Hie one dimensional case. One 
can expect a scale free behaviour when the network shows 
small world effect. The absence of scale free behaviour 
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